Van-der-Waals supercritical fluid: Exact formulas for special lines 
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In the framework of the van-der-Waals model, analytical expressions for the locus of extrema 
(ridges) for heat capacity, thermal expansion coefficient, compressibility, density fluctuation, and 
sound velocity in the supercritical region have been obtained. It was found that the ridges for 
different thermodynamic values virtually merge into single Widom line only at T < 1.07T C ,P < 
1.25P C and become smeared at T < 2T C , P < 5P C , where T c and P c are the critical temperature and 
pressure. The behavior of the Batschinski lines and the pseudo-Gruneisen parameter 7 of a van-der- 
Waals fluid were analyzed. In the critical point, the van-der-Waals fluid has 7 = 8/3, corresponding 
to a soft sphere particle system with exponent n = 14. 

PACS numbers: 64.10. +h, 64.60.F-, 65.20.De 



A liquid-gas phase equilibrium curve onto the T,P - 
plane ends at the critical point. At pressures and tem- 
peratures above critical ones (P > P C ,T > T c ), the prop- 
erties of a substance in the isotherms and isobars vary 
continuously, and it is commonly said that the substance 
is in its supercritical fluid state, where there is no dif- 
ference between a liquid and a gas. An anomalous be- 
havior of the majority of characteristics are observed in 
the vicinity of the critical point. The correlation length 
for thermodynamic fluctuations diverges at the critical 
point one can also observe a critical behavior of the 
compressibility coefficient /3y, thermal expansion coef- 
ficient ap, and heat capacity C p : the given properties 
pass through their maxima under a change of pressure or 
temperature. Near the critical point, the positions of the 
maxima of these values in the T, P - plane are close to 
each other [![. The same is true for the value of density 
fluctuations, the speed of sound, thermal conductivity, 
etc. Thus, in the supercritical region, there is a whole 
set of the lines of extrema of various thermodynamic val- 
ues. Each of these lines can be regarded as a "thermo- 
dynamical" continuation of the liquid-gas phase equilib- 
rium curve into the supercritical region. The smearing 
and decreasing (in magnitude) extrema of each of the val- 
ues form a " ridge" [2t4| . A knowledge of the positions of 
the above "ridges" in the T,P - plane is very important; 
in particular, it determines a maximum value for such 
technologically essential characteristics as the dissolving 
ability of a supercritical fluid, the rate of chemical reac- 
tions in a fluid, and others (0^3 and refs therein). It 
turned out that the experimentally observed lines of the 
"ridges" are close to an isochores with a slight decrease 
in density with increasing temperature 0-0] . Most stud- 
ies on the supercritical region focused on examining the 
"ridge" for the density fluctuations 

G. Stanley suggested the name "Widom line" for the 
line of the maximum of the correlation length isotherms 
and isobars Since the lines of the maxima near the 
critical point merge into one line, the above term was 



proposed to be used in a wider sense - in reference to 
the lines of the maxima of all values determined by the 
second derivatives of the Gibbs thermodynamic poten- 
tial. The Widom lines for liquid - liquid and more rarely 
for liquid - gas transitions for a number of systems were 
considered (5M9(. It is a priori unclear how far from the 
critical point we may say of a single Widom line for all 
ridges, and how far from the critical point the extrema of 
particular physical quantities can still be followed. The 
first part of question may be rephrased as " how far from 
the critical point the prefactor that appears in the ex- 
pression of fluctuations, and response functions, in terms 
correlation length is close to a constant" . 

Apart from the Widom line, other "special" lines, sep- 
arating a fluid state, have been suggested. The Batschin- 
ski line [1(1, sometimes referred to as the Zeno line 



111 ll 21 ] , corresponds to a formal coincidence between the 
equation of state for a fluid and the equation of state for 
an ideal gas. 

For solids, an important thermodynamic parameter is 
the Gruneisen parameter which reflects the variation of 
the lattice dynamics under a change of density. For a 
fluid, one can introduce a pseudo-Gruneisen parameter 
[l3j , relating such thermodynamic quantities as the heat 
capacity C v , thermal expansion coefficient ap and com- 
pressibility coefficient Pt- 



y = a P /(C v p T ). 



(1) 



It is of great interest to analyze the behavior of "spe- 
cial" lines in a fluid and the pseudo-Gruneisen parameter 
in the framework of a simple, exactly solvable model. In 
the present study we have analyzed the properties of a 
van-der-Waals fluid. The van-der-Waals equation is one 
of the simplest equations of state for a fluid. In the re- 
duced variables T r — T/T c ,P r = P/P c ,p r = p/p c the 
equation has the form: 



(P r + 3p 2 r )(3- Pr ) = 8T r p r . 



(2) 



It is well known that the van-der-Waals equation does 
not reproduces the behavior of real fluids completely and 
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FIG. 1: (Color online) Positions of the maxima of thermal 
expansion coefficient ap, compressibility Pt, heat capacity 
Cp and density fluctuations C,t in the (p — T) (a) and (P — T) 
(b) coordinates. The arrow indicates the approximate end of 
the "single" Widom line. Thick lines correspond to liquid-gas 
transition. 



precisely ( [1 , 14 1 and refs therein) . At the same time the 
van-der-Waals equation can be used to understand the 
fundamentals of fluid behavior. Besides, this equation is 
advantageous because exact analytical expressions can be 
obtained for most of the physical quantities [l4j . At the 
same time, however strange it may look, a supercritical 
region of the T,P - parameters for the van-der-Waals 
fluid model has been studied insufficiently. We know of 
only one study [H analyzing the behavior of the line of 
the maximum of density fluctuations < (AN) 2 > / < 
N >= T(dp/dP) T = Ct- It is found that the position 
of the line of the maxima of density fluctuations on the 
isotherms satisfies the equation: 



P, 



= 3-271 



1/3 



(3) 



Thus, the above line formally ends at zero density and 
zero pressure at T r — 27/8 (see Fig. [T]). The speed of 
sound has the form: 



v s = ^/t/Ct = 



/ 6(p r (3-p r ) 2 -4r r ) 

\ p r (3 - p r ) 2 



1/2 



The line of the minimum of the speed of sound on the 
isotherms obviously corresponds to the equation (J3j) . In 
it was supposed that the extrema of other thermo- 
dynamic quantities in the isotherms would lie roughly 
on the same line. However, as we will show below, all 
" ridges" diverge as one goes even slightly away from the 
critical point. 

Isothermal compressibility in the framework of the van- 
der-Waals model has the form: 



dp 

dP 



(Pr ~ 3) 2 



6pr(— 4T r + p r (p r ~ 3) 2 ) 



The line of the maxima of compressibility Pt on the 
isotherms satisfies the equation: 



T r = 



p r {3 - p r ) 3 

2(3 + Pr ) 



(4) 



This line ends at its own critical point at T r — 1.156, P r — 
1.285, p r = 0.646 (see Fig.Q]). Using a known thermody- 
namic relation [![, we obtain: 



Cp — Cy 



32T r 



3(4T r - 9p r + 6p 2 - pf) 



We remind that at T r > 1 CV = 3/2. It can easily be 
seen that the line of the maxima of the heat capacity Cp 
in the isotherms coincides with the isochore p r = 1 and 
is described in the T, P coordinates by the equation: 



_ 3 1 

T r - - + -P r 



(5) 



i.e., represents a direct continuation of the gas-liquid 
equilibrium line T| (see Fig. [T]). The thermal expan- 
sion coefficient has the form: 



ap 



dp 
df 



4(p r - 3) 



3p r {p r -3) 2 - 12T r 



The line of the maxima of the thermal expansion coeffi- 
cient ap in the isotherms corresponds to the equation 



Tr 



(3-2p)(3-p) 2 /4. 



(6) 



This above line ends at zero pressure and zero density at 
T r = 27/4 (see Fig. [p. 

Although all ridges are described by different equa- 
tions, they are close together near the critical point. For 
the estimate, the lines of the extrema can be thought as 
coinciding, if the temperature values on the lines at the 
same pressure differ by less than 1%. The value of 1% 
roughly corresponds to the experimental accuracy of a 
measurement of the respective values and to the errors in 
the computer simulation data. For van-der-Waals fluid 
the positions of all ridges for different thermodynamic 
values merge into " single" Widom line in the P, T- co- 
ordinates at T < 1.07T C , P < 1.25P C (see the inset in 
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Fig. O^b)). Although the line of the maxima of density 
fluctuations and that of the maxima of thermal expan- 
sion coefficient for a van-der-Waals fluid formally go till 
zero pressures, and the line of the maxima of heat capac- 
ity goes to an infinite temperature region, the amplitude 
of all extrema very rapidly decays if we move away from 
the critical point (see Fig. [2]). It is of interest to notice 
the existence of an additional minimum on the tempera- 
ture dependence of compressibility 0t in the isobars (see 
the inset in Fig.^b)). 

As a criterion of actual disappearance of the extremum, 
one can consider the ratio of a respective thermodynamic 
value in the maximum or minimum to this value at den- 
sities being 10% different from the density in the ex- 
tremum. If this ratio is below 1.01 (the difference be- 
tween the extremal value and the "background" value 
is below 1%), the "ridge" can be thought of as actually 
smeared. When using the above criterion, the lines of all 
extrema, in fact, end at rather moderate temperatures 
and pressures: T r ~ 1.59, P r ~ 2.78, p r ~ 0.83 for ap; 
T r ~ 1.44, P r ~ 2.13, p r ~ 0.74 for ( T ; T r ~ 1.73, P r ~ 
3.9, p r = 1 for C p ; T r ~ 1.15, P r ~ 1.35, p r ~ 0.73 for 
P T ; T r ~ 1.28, P r ~ 1.84, p r ~ 0.83 for V s (see inset in 
Fig. ttTb)). 

The lines of the maxima of most quantities correspond 
to a decrease in density with increased temperature (see 
Fig. [1]); only the line of the maxima of the heat capacity 
lies on the isochore. 



FIG. 3: (Color online) The Batschinski lines in the (p-T) (a) 
and (P — T) (b) planes. Thick lines correspond to a liquid-gas 
transition. 

Thus, a "thermodynamic" continuation of the gas- 
liquid phase equilibrium line for a van-der-Waals fluid 
represents single Widom line if the temperature moves 
only 7 % away from the critical point; if the tempera- 
ture moves further away, it represents a rapidly widening 
bunch of lines, which in fact ends at T r < 2 P r < 5. 

Let us discuss now the Batschinski line, sometimes re- 
ferred to as the Zeno line. This line corresponds to the 
equation PV/RT = 1. A.I. Batschinski demonstrated 
that within the van-der-Waals model, the above line will 
be straight in the coordinates p, T [10(. In most studies, 
the behavior of this line was only analyzed in the p, T 
coordinates for model and real systems [ill, E3] ■ At the 
same time it is of interest to examine not only the be- 
havior of this line but also the behavior of other lines 
determined by the condition PV/RT = Z both in the 
p, T and T, P-planes. Fig. [3] presents the lines satisfying 
the conditions PV/RT — Z for a van-der-Waals fluid, 
obtained from the equation: 

_ 3p r (3 - p r ) 
r ~8-k(3-p r )- [) 

For the Batschinski line, we have the expression T r = 
| (3 — p r ). In this case, it should be taken into ac- 
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count that the condition PV/RT = Z in the reduced 
units has the form: ~pr~ — where k = 8/3 corre- 
sponds to the condition PV/RT = 1. The behavior of 
the lines if PV/RT < 1 and PV/RT > 1 is quite differ- 
ent. The Batschinski line, determined by the equation 
PV/RT = 1, separates two regions of a fluid: a "soft" 
low density fluid with the predominance of the attractive 
potential as compared to an ideal gas, where PV/RT <1, 
and a more "rigid" higher density fluid with the predom- 
inance of the repulsive potential as compared to an ideal 
gas, where PV/RT > 1 (see Fig. [3]). Therefore, the line 
PV — RT can be called " separatrix" . The Batschin- 
ski separatrix is the only line of the given family which 
ends at zero density and pressure, where it coincides with 
the ideal gas regime. The Batschinski separatrix for the 
van-der-Waals fluid ends at T r = 27/8, i.e., the same 
temperature at which the lines of the maxima of density 
fluctuations and minima of the speed of sound end. This 
coincidence is related to the fact that these lines corre- 
spond to the zero of the second pressure derivative of the 
density. For an ideal gas, this derivative is equal to zero 
at all temperatures, so the line of the maxima of the value 
at zero density corresponds to the ideal gas equation 
PV = RT as well. 

Finally, let us analyze the behavior of a pseudo- 
Gruneisen parameter for a van-der-Waals fluid. Using 
the equation ^ we have: 
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16 p r 
3 (3 - p r ) ' 



(8) 



For the majority of real substances, the Gruneisen pa- 
rameter varies in the range from 0.5 to 3, i.e., the van-der- 
Waals fluid in the region of "normal" densities (2-2.5) for 
liquids has anomalously high Gruneisen parameter val- 
ues. For a system of particles with a purely repulsive ex- 
ponential interaction (soft sphere system), the Gruneisen 
parameter value can easily be deduced: 7 = (n + 2)/6, 
where n - is the exponent of the repulsive potential. 
Thus, p r = 3 limit obviously corresponds to hard sphere 
system (n = 00). At the critical point, 7 = 8/3 « 2.67, 
that is, coincides with the value for a for soft sphere sys- 
tem with n — 14, which is close to the exponent in the 
repulsive part of Lennard-Jones potential (n = 12). The 
behavior of real rare gas substances, too, is well described 
by the potential of soft spheres with n — 12 — 13. Thus, 
the melting curve for argon coincides, to a high accu- 
racy, with that for a soft sphere system with n = 12.7 
[lif . Rare gas solids and fluids also have typ ical values 
of the Gruneisen parameter 7 ~ 2.3 — 2.6 Jjjlla]. Thus, 
the van-der-Waals fluid near the critical point has the 
Gruneisen parameter values close to those for real rare 
gas fluids. In this regard, a success in the description of 
the properties of rare gas fluids by van-der-Waals equa- 
tion in the vicinity of the critical point can better be 
understood. 



Summing up, one can conclude that in the framework 
of the van-der-Waals equation it has been possible to ob- 
tain exact analytical expressions for the "special" lines 
in the region of a supercritical fluid and for the pseudo- 
Gruneisen parameter. The qualitative inferences made in 
the present study will obviously be valid for real simple 
fluids as well. This is particularly true of the conclusion 
about how far from the critical point we may establish 
a single Widom line for different thermodynamic values, 
and how far from the critical point the extrema of partic- 
ular physical quantities can still be followed. Of course, 
the behavior of ridges for thermodynamic values may be 
different and depends on corresponding equation of state 
However, as for the liquid-gas transitions at enough 
high pressures (P r > 10), it is dynamic, not thermody- 
namic, characteristics that should be considered when 
separating a fluid into liquid-like and gas-like regions. 
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